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ABSTRACT 

We present a physical model for the origin of z ~ 2 Dust-Obscured Galaxies (DOGs), 
a class of high-redshift ULIRGs selected at 24 |j.m which are particularly optically faint 
(-p24jum/ Fr > 1000). By combining iV-body/SPH simulations of high redshift galaxy evo- 
lution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs 
(with F-2A> 0.3 mJy at z ~ 2) are well-modeled as extreme gas-rich mergers in massive 
(~ 5 x 10 12 - 10 13 Mq) halos, with elevated star formation rates (~ 500 - 1000 MQyr" 1 ) 

and/or significant AGN growth (Mbh ^ 0.5 MQyr -1 ), whereas less luminous DOGs are 
more diverse in nature. At final coalescence, merger-driven DOGs transition from being star- 
burst dominated to AGN dominated, evolving from a "bump" to a power-law shaped mid-IR 
(IRAC) spectral energy distribution (SED). After the DOG phase, the galaxy settles back to 
exhibiting a "bump" SED with bluer colors and lower star formation rates. While canonically 
power-law galaxies are associated with being AGN-dominated, we find that the power-law 
mid-IR SED can owe both to direct AGN contribution, as well as to a heavily dust obscured 
stellar bump at times that the galaxy is starburst dominated. Thus power-law galaxies can be 
either starburst or AGN dominated. Less luminous DOGs can be well-represented either by 
mergers, or by massive (Afbaryon ~ 5 x 10 11 Mq ) secularly evolving gas-rich disc galaxies 
(with SFR > 50 MQyr -1 ). By utilising similar models as those employed in the SMG for- 
mation study of Narayanan et al. (2010), we investigate the connection between DOGs and 
SMGs. We find that the most heavily star-forming merger driven DOGs can be selected as 
Submillimetre Galaxies (SMGs), while both merger-driven and secularly evolving DOGs typi- 
cally satisfy the BzK selection criteria. The model SEDs from the simulated galaxies match ob- 
served data reasonably well, though Mrk 231 and Arp 220 templates provide worse matches. 
Our models provide testable predictions of the physical masses, dust temperatures, CO line 
widths and location on the MBH-Afbuige relation of DOGs. Finally, we provide public SED 
templates derived from these simulations. 

Key words: cosmology:theory-galaxies:formation-galaxies:high-redshift- 
galaxies:starburst-galaxies:ISM-galaxies:ISM-ISM:dust 
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1 INTRODUCTION 

Redshift z ~ 2 is a rich epoch for understanding galaxy forma- 
tion and evolution. During th is time period, the bulk of the cosmic 
stellar mass was assembled (Dickinson et alj|2003t I Rudni ck et al.l 
2006), and the star forma tion rate and quasar space densities 
are both near their peaks (|Bouwens et al.l 120041 ; IHopkins! |2004| ; 
IHopkins et alj|2007l : lshaver et aljl996h~ 

In the local Universe, infrared (IR) selection of galaxies is an 
efficient means of selecting galaxies undergoing significant star for- 
mation (e.g., Arp 220), and possibly hosting optically-obscured ac- 
tive galactic nuclei (AGN; e.g., Mrk 231, NGC 6240). By anal- 
ogy, at higher redshifts, IR selection techniques select galaxies 
which are major contributors to the cosmic star formation rate 
density, far-IR background, and progenitors of the present-day 
massive galaxy population. Indeed, ultraluminous infrared galax- 
ies (ULIRGs; Lir> 10 12 Lq) at z ~ 2 app ear to contribute sub- 
stantially to the infrared l uminosity density jPerez-Gonzalez et al.l 
2001 ICaputi et alj|2007l; iReddv et all 120081 : IHopkins et alJtoTa 
Hopkins & Hernquist 2010), and the z>0.7 star formation rate den- 
sity dLe Floc'h et alj|2005l ; lshi et alj2009l) . 

A well-studied sample of galaxies selected for their FIR prop- 
erties at z ~ 2 are Submillimetre Galaxies (SMGs), chosen in 
blank field surveys for their redshifted cool dust emission at 850 
(j.m above F$5o > 5 mly. These galaxies are in a transient star- 
burst phase (e.gJSwinbank et alj 2004 ; Menende z-Delmestre et al.l 



120071: 1 Younger et alj|2008 j ; iMichatowski et alj|2009l). host rapidly 
growi ng supermassive black holes (e.g.. lAlexander et alj 120051 
2008) , and resid e in extremely massive -~5x 10 12 M halos 



i Blain et al.l |2004> . While these galaxies represent an important 
sub-population of infrared-luminous galaxies at the epoch of 
peak galaxy formation, there is some concern that a selection at 
850 (im may be missing significant populations of high-redshift 
ULIRGs with warmer dust temperatures (e.g., iBussmann et al 



2009 



2009 



ist tem pe 
]2009t |C 



Casey et alj 12009k iHuang et al 



jChapman & Casevl p 

lYounger et alj 12009). As such, alternative selection tech- 
niques for identifying z ~ 2 ULIRGs have become desirable. 

With the launch of the Spitzer Space Telescope, surveys at 
24 yim have uncovered a population of ULIRGs at z ~ 2 (e.g., 



Rigbv et al. 200' 
2008; Soiferetal 



Donley et alj2007l;lYan et al 



2008; lLonsdale et ail 2009: 



2007; Farrah et al, 



Huang et alj|2009l 



and refe rences therein). Re cently jDev et ak d2008[ hereafter. D08 
see also lFiore et al. I d2008l) ) presented an efficient means of identi- 
fying high-redshift ULIRGs which are both mid-infrared luminous 
and optically faint. Specifically, by imposing a nominal selection 
criteria of i?-[24] > 14 (roughly corresponding to a flux density 
ratio F 24 /Fr> 1000; formally F 24 /F R > 960) and F 24 > 0.3 mjy, 
D08 found a sample of ULIRGs with a relatively narrow redshift 
distribution (centered at <z>=1.99 with ct z =0.5) which exhibit 
heavy reddeni ng of the intrinsic UV flux and str ong rest-frame 8 
Urn emission dHouck et ai]|2005l ; iFiore et alJbooj) . D08 refer to 
these galaxies as Dust-Obscured Galaxies (DOGsQ. 

DOGs exhibit a range of luminosities, though AGN- 
dominated galaxies tend to preferentially d ominate the bright 



end of the 24 pm flu x density distribution (B rand et al 
Weedman etai]|2006j|bl ; iDev et alj|2008l ; IFiore et alj|200~ 



2006 



2009 



Sacchi et alj|2009r) . The DOGs with the highest 24 p.m flux den- 



1 We note that while this paper nominally focuses on 24 pm-selected 
galaxies which additionally have 24/i?> 1000, the results of this paper are 
generally applicable to the 24 pm-bright (-F24> 0.3mJy) z ~ 2 ULIRG 
population. 



sities (f24> 0-8 mjy) have bolome t ric infrared lumin osities (Lir) 
~ 10 13 Lq teussmann ei"al]|2009l ; iTvler et alj|2009l) . As a popu- 
lation, the subset of F 2 4> 0.3 mly galaxies which are also DOGs 
contribut e around a quarte r of the total infrared luminosity density 
at z ~ 2 jPev et alj|2008l) . and constitute ~60% of the total con- 
tribution by z ~ 2 ULIRGs. It is clear that DOGs are a cosmolog- 
ically significant population of galaxies with diverse observational 
characteristics. 

Broadly, 24 pm-selected ULIRGs at z ~ 2 (including DOGs) 
appear to fall into two categories based on their mid-IR IRAC 
SEDs: those with a power-law shape and those which exhibit a 
bump at observed A ~ 5 pm corresponding to the 1.6 pm rest 
frame stellar photospheric bump (hereafter, these are referred to as 
PL galaxies and bump galaxies, respectively). Bump galaxies are 
thought to be dominated by star formation, with the bump (at rest- 
frame 1.6 [im) originating from starlight. PL galaxies are thought 
to ha ve their mid-IR SEP dominated by AGN continuum emis- 
sion dWeedman et al.l l2006bl ; IPonlev et al]|2007l: lYan et alj|2007l; 
iMurphv et alj|2009h . In support of this scenario, PL galaxies have 
SEPs whi ch are reasonably represented by that of a Mrk 231 
template dBussmann et alj 120091 ; ITvler et alj 120091) . and tend to 
show relatively compact opt i cal morphologies (R e ~ 1 — 6 kpc; 
Melbourne etai]|2008L I2OO9I ; IBussmann et al.l l2009). On the other 
hand, bump galaxies are typically polycyclic aromatic hydrocar- 
bon (PAH)-rich, show PAH e quivalent widths consistent with be- 
ing star formation dominated dYan et~a l. 2005; Saiin a^et"aT]|2007l ; 
lFarrahetalJl2008l ; IPesai et alj|2009l: IHuang et alJbOOgTT and are 
more often morpholog ically resol ved into multiple compo nents 
( IBussmann et alj [2009f) . However, IFiore et ak I j2008L l2009h have 
argued that even these less luminous galaxies may harbor heavily 
obscured AGN. 

DOGs are massive galaxie s , with stellar masses M* > 
10 10 " 11 M q (IBussmann et al.l 120091 ; lLonsdale et all 120091; 
IHuang et alj 120091) . and cluster in group-siz ed halos of order 
(M DM » 10 12 " 13 Mq ; iBrodwin et alj[20o3) . The clustering is 
luminosity-dependent, with more luminous DOGs inhabiting more 
massive halos. Since the more luminous DOGs also tend to exhibit 
power-law SEDs, it follows that PL DOGs cluster more strongly 
than the lower luminosity bump DOGs. 

Although much progress has been made in understanding 
DOGs in a relatively short time period, a myriad of fundamental 
questions regarding their physical nature exist. Are DOGs prefer- 
entially isolated galaxies or mergers (or both)? Does the 24/ R cri- 
terion select galaxies at a particular evolutionary point? How are 
bump DOGs and PL DOGs related - are they distinct galaxy popu- 
lations, or possibly related via an evolutionary sequence? How are 
DOGs related to other coeval high-redshift populations (e.g., the 
SMGs, BzK galaxies, etc.)? 

Numerical simulations can offer complementary information 
to the observations in hand, and provide insight into the afore- 
mentioned questions. In particular, by coupling radiative transfer 
modeling with hydrodynamic simulations of galaxy evolution, one 
can precisely mimic local and high-redshift observational selec- 
tion functions, and directly relate obs erved trends to physical con- 



2009) 



ditions in the model galaxies (e.g., Ipnsson 2006; lonss on et al 



20061 l20ld: lLotz et aU2008l : iNaravanan et al 1120 id : lYounger et al 



Previously, we have utilised similar methods to investigate the 
formation and evolution of h igh-redshift Submillimetre Galaxies 
INaravanan eTai] d2009L I20T0I C. Hayward et al. in prep.). Here, 
we build up on these efforts by util ising similar models as thos e 
employed in Narayan an et alj d2009l) and INaravanan et all d2010l) . 
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This will allow us not only to investigate a potential formation 
mechanism for DOGs, but to explore the potential connection be- 
tween z ~ 2 ULIRGs selected for their 24 |j.m brightness, and those 
selected in the submillimetre. 

The paper is organised as follows. In §[2] we detail our numer- 
ical methods. In § [3] we describe the formation of z ~ 2 DOGs, 
and explain the evolution of the 24/R ratio and the relationship be- 
tween bump and PL DOGS. In §|4] we discuss the physical require- 
ments necessary to form a DOG, and whether mergers are neces- 
sary; In § [5] we detail the relationship between DOGs, SMGs and 
BzK galaxies. In §[6] we relate our model to existing observations, 
and make testable predictions. In §[7] we provide discussion. We 
conclude in §[8] 

Throughout this paper, we utilise a fiducial DOG selection cri- 
teria 24/R > 1000 and F 24 > 300 /iJy, consistent with the DOGs 
sample in D08. We note again that while we focus on these 24 
\im selected galaxies which are optically faint (DOGs), the results 
are generally applicable to most 24 |j.m-selected z ~ 2 ULIRGs. 
Throughout, we assume a cosmology with Q,a = 0.7, Q, m — 0.3 
and h = 0.7. Finally, we note that the models presented here are not 
cosmological in nature. As such, all non-evolutionary plots should 
be taken to represent ranges of expected values and colors, rather 
than true expected distributions. 



2 NUMERICAL METHODS 

Our methodology is as follows. We first run a suite of Smoothed- 
particle hydrodynamics (SPH) simulations of binary gas-rich 
galaxy mergers over a range of baryonic masses and merger mass 
ratios. The mergers are between idealised disc galaxies. We ad- 
ditionally run three isolated disc galaxy simulations as a control. 
These simulations represent potential physical models for z ~ 2 
24 |j.m sources. We then investigate the UV-mm wave SEDs of 
these galaxies utilising a 3D dust radiative transfer code in post- 
processing. In what follows, we describe the details of these nu- 
merical simulations, and outline both the input parameters, as well 
as the assumptions made. 



2.1 Hydrodynamics 

The hydrodynamic evolution of the galaxies was simulated utilis- 
ing a modifi ed version of th e publicly available iV-body/SPH code, 
GADGET- 2 ( Springe] 2005). GADGET- 2 follows the dynamical evo- 
lution of the dark matter, ISM, stars and black holes in the simu- 
lated galaxies, and explicitl y conserves both energy and entropy 
dSpringel & Hernquistll2002h . 

The ISM is c onsidered to comprise of a tw o-phase medium via 
a sub-grid model (Springel & Hernquist 2003). In this model, cold 
clouds are assumed to be embedded in a hotter, pressure-confining 
phase. The cold phase may convert to the hotter phase via evapora- 
tion from supernovae heating, and similarly the hot phase may cool 
into cold gas. This ISM is pressurised by heating from supernovae. 
Numerically, this is handle d via an effective equa tion of state (here- 
after EOS; see Figure 4 of ISpringeletall20 05b). Here, we assume 
the stiffest equation of state which stabilises the disc against run- 
away fragmentation. We explore the consequences of relaxing this 
assumption in § 12.2.21 

The star formation rate in the simulations is dependent on the 
density of cold gas (defined here as gas which has reached the floor 
of the cooling curve, T — 10 4 K) as SFRoc Pg a g. The normalisation 
of this relation is set such that the simulated star formation rate 



surface density scales with the gas surface density in agreement 
with local observations dKennicuttlll998l ; ICox et alJuOOd) . 

Black holes are included in the simulations as collisionless 
particles. These particles stochastically accrete matter according to 
a Bondi-Hoyle-Lyttleton paramaterisation when a particle comes 
within a smoothing length of the black hole sink particle. To model 
the impact of thermal feedback from accreting black holes, a frac- 
tion of the accreted mass energy (here, 0.5%) is r einjected ther- 
mally and spherically into the surrounding ISM ( Sprin ge! et al.l 
|2005bj). The fraction of accreted mass energy which couples back 
to the ISM is set by matching the normalisation to the local 
Mbh — o relation (though n ote that the slope is a natural conse- 
quence of the AGN feedback; |Pi Matteo et alj2005h . 

The progenitor disc galaxies are initialised with a [Hernquist 

dl990h dark matter halo profile. The choice of concentration and 
virial radius for a halo of a given mass was motivated by cos- 
mological TV-bod y simulations, then scaled to match the ex pected 
redshift evolution dBullock et al.ll200ll ; lRobertson et aT]|2006b . The 
halo for our fiducial model (DOG10, Table[TJ has a z=0 concentra- 
tion c = 11.0 and spin parameter A = 0.033. We initialise our discs 
to reside at z=i, so that our simulated mergers occur around z ~ 2. 
We consider initial virial velocities in the range 225-500 km s , 
resulting in dark matter halo masses of ~ 10 12 Mq- IO^Mq, con- 
sistent with the DOG clustering measurements by iBrodwin et al.l 
(2008). The halos are popula ted with ex ponential discs whose disc 
scale lengths follow the lMo et alj dl998h prescriptions. 

The discs for model galaxies are initialised with gas frac- 
tions / g = gas /(ga s + stars) = 0.8. Simulation tests by 



INaravanan et alj ( 2010) have shown that high-redshift mergers ini- 
tialised with this gas fraction result in / g ~ 0.2 — 0.4 at the 
time when the galaxy may be viewed as an SMG, which is com- 
parable to observations (e.g Bouche et alj 1 20071 ; iTacconi et alj 
120081: iDaddi et alj|2009l : ITacconi et alj|2010h . While gas fractions 
in DOGs are observationally unconstrained, based on an observed 
overlap between the luminous DOG and SMG populations (e.g., 
|Popeetal.ll2008l ;l Brod win et "al] |2008l) . we consider this assump- 
tion reasonable. In order to investigate any potential dependence 
on gas fraction, we have run a subset of our simulations with initial 
gas fractions f g — 0.4, 0.6 as well. The results from all initiali- 
sations are qualitatively similar, and we include both in our plots 
to increase the simulation sample size. The gravitational soften- 
ing lengths are set at 100 h _1 pc for baryons and 200 h _1 pc for 
dark matter. Because structure in the ISM clearly exists below these 
scale lengths, we employ a subgrid specification to the ISM for the 
radiative transfer calculations. This will be discussed in the follow- 
ing two subsections. 

In this paper, we consider 16 simulations, with inital con- 
ditions listed in Table JT] The simulatio ns (and nomenc lature) 
here are nearly identic al 2 ] to those in INaravanan et alj (120090 and 
INaravanan et alj d2010h 's studies of SMG formation. This will fa- 
cilitate direct comparisons between the DOG population and high- 
redshift SMGs. When we refer to investigating a fiducial merger, 
this will refer to model DOG10 with initial f g = 0.8. 



2 The pr imary difference betwe en the simulations run here and those em- 
ployed in INaravanan e7aH (2009) is in the assumptions of the evolving dust 
mass (see the next section). Here, we tie the dust mass to the metallicity, 
whereas in Narayanan et al ] <2009l) . the dust mass is tied to the gas mass. 
This makes little difference during the merger event, though can have some 
effects on the optical flux during the inspiral phase. 
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Table 1. DOGs are ordered with decreasing halo mass (and grouped by mergers and 
isolated galaxies, such that the isolated galaxies have the prefix 'i' in their name). Column 
1 is the name of the model used in this work. We emphasise that these model s are n early 
identical to those employed in the SMG formation study by Naray anan et alj <2010h and 
iNaravan an et al. 1 2003). This will facilitate the comparison of these models with SMGs 
(ij |5. 1 \ . Columns 2 and 3 give the virial velocity and halo mass of the galaxies. Column 4 
is the total baryonic mass of the system. Column 5 is the merger mass ratio. Columns 6 & 
7 are initial orientations for disc 1 , Columns 8 & 9 are for disc 2. The orientation angles 
are with respect to the merger plane (such that 9,<f> = would be a coplanar merger). 
Column 10 lists the initial gas fraction of the simulation. 



Model 


V c 
(kms -1 ) 


A/dm 
M 


A/ bar 
M Q 


Mass Ratio 


01 


01 


82 


4>2 


f 3 


DOG1 


500:500 


1.25xl0 13 :1.25xl0 13 


l.lxlO 12 


1:1 


30 


60 


-30 


45 


0.8,0.6,0.4 


DOG2 


500:500 


1.25xl0 13 :1.25xl0 13 


1.1 xlO 12 


1:1 


360 


60 


150 





0.8 


DOG3 


500:500 


1.25xl0 13 :1.25xl0 13 


1.1 xlO 12 


1:1 


-109 


-30 


71 


-30 


0.8 


DOG4 


500:320 


1.25xl0 13 :3.4xl0 12 


6.9X10 11 


1:3 


30 


60 


-30 


45 


0.8,0.6,0.4 


DOG5 


500:320 


1.25xl0 13 :3.4xl0 12 


6.9xlO n 


1:3 


360 


60 


150 





0.8 


DOG6 


500:320 


1.25xl0 13 :3.4xl0 12 


6.9X10 11 


1:3 


-109 


-30 


71 


-30 


0.8 


DOG7 


500:225 


1.25xl0 13 :1.2xl0 12 


6.0X10 11 


1:10 


30 


60 


-30 


45 


0.8,0.6,0.4 


DOG8 


500:225 


1.25xl0 13 :1.2xl0 12 


6.0X10 11 


1:10 


360 


60 


150 





0.8 


DOG9 


500:225 


1.25xl0 13 :1.2xl0 12 


6.0X10 11 


1:10 


-109 


-30 


71 


-30 


0.8 


DOG 10 


320:320 


3.4xl0 12 :3.4xl() 12 


2.9 xlO 11 


1:1 


30 


60 


-30 


45 


0.8,0.6,0.4 


DOG 11 


320:320 


3.4xl0 12 :3.4xl0 12 


2.9X10 11 


1:1 


360 


60 


150 





0.8 


DOG12 


320:320 


3.4xl0 12 :3.4xl0 12 


2.9xlO n 


1:1 


-109 


-30 


71 


-30 


0.8 


DOG 13 


225:225 


1.2xl0 12 :1.2xl() 12 


l.OxlO 11 


1:1 


30 


60 


-30 


45 


0.8,0.6,0.4 


iDOGl 


500 


1.25X10 13 


5.5 xlO 11 




N/A 


N/A 


N/A 


N/A 


0.8,0.6,0.4 


iDOG2 


320 


3.4xl0 12 


1.5xlO n 




N/A 


N/A 


N/A 


N/A 


0.8,0.6,0.4 


iDOG3 


225 


1.2xl0 12 


5xl0 10 




N/A 


N/A 


N/A 


N/A 


0.8,0.6,0.4 



2.2 Radiative Transfer 

2.2.1 General Methods 

We extract synthetic SEDs from our SPH simulations in post- 
processing utilising the 3D polychromatic Monte Carlo radia - 
tive transfer code, SUNRISE djonssoni 2006; Jon sson et aDl2010h . 
SUNRISE considers the propagation of UV through mm-wave pho- 
tons through a dusty medium, and operates on an adaptive mesh. 
Here, we summarise the as pects of the code mo st relevant to this 
study, and refer the reader to Jon sson et al. for the most cur- 

rent detailed description of the underlying algorithm. 

Dust radiative transfer may be thought of in terms of sources 
and sinks. Here, the sources of radiation are stellar clusters, AGN 
and dust grains. The stellar clusters emit a spectrum calculated b y 
STARBURST99 iLeifherer et alJll999l : IVazquez & Lei fherel 2005). 
where the ages and metallicities of the stars are known from the 
SPH simulations. The stellar particles present in the simulation's 
initial conditions are assumed to have formed over a constant star 
formation history for ~250 Myr (where this time was chosen by di- 
viding the stellar mass of the first simulation snapshot by the SFR 
at that time). That said, the rest-frame UV and mid-IR (which dom- 
inates the wavelength regime of interest for DOGs) is contributed 
mainly by young (< 100 Myr) stars (and their surrounding HII re- 
gions), and thus our results are largely insensitive to the choice of 
star formation history. We assume a Rroupa IMF f or these calcu- 
lation s, suggested by observations of z ~ 2 SMGs jTacconi et al.l 
l2008h . 

The A GN input spectrum u tilises the luminosity-dependent 
templates o f iHopkins et afl J20071) of type I quasars. The template 
spectrum includes to rus hot dust emi s sion, a nd is based on the mid- 
IR SED template of iRichards et al.l J2006l) . The normalisation of 
the input spectrum is set by the total bolometric luminosity of the 



black hole(s), which is determined by their mass accretion rate: 
£agn = ijMbhc 2 . Here, r\ is assumed to be 10%. The AGN are 
included in all simulations investigated here (mergers as well as 
secularly evolving galaxies). 

The sinks of radiation are dust grains in the ISM. The emission 
from stars and AGN undergo scattering, absorption and remission 
as they propagate through the ISM. Because the ISM substructure is 
not resolved by our hydrodynamic simulations, we employ subgrid 
modeling techniques which we explore in the following subsection. 
Eight model cameras are set up isotropically around the galaxy to 
observe the merger. The emergent flux is determined by the number 
of photons that escape the galaxy in a given camera's direction. 

SUNRISE calculates the dust temperature of the grains self- 
consistently, and assumes that the dust and radiati on field are in 
radiative equilibrium utilising the methodology of ljuvelal J2005h 
for calculating the converged radiation field. In this methodology, 
when photons are absorbed in a grid cell, the dust temperature is 
updated. A new photon packet is then emitted from this cell which 
has an SED equal to the difference between the SED associated 
with the updated dust temperature, and that from the old one. This 
procedure is iterated upon until the radiation field has converged 
(Jonsson et al. 2010). Consequently, both the AGN and the star- 
burst are involved in determining the dust temperature across all 
the grid cells. 

PAHs are emitted by the HII regions, PDRs, and the diffuse 
dust in the galaxy. Carbonaceous dust grains with sizes a < 100 A 
are assumed to have characteristics of PAH molecules. A fraction 
of these PAH molecules (following Ijonsson et alJl201cL we take 
this value to be 50%) emit a template spectrum, and the remain- 
ing fraction emit thermally. The PAH template is a linear combina- 
tion of Lorentzian profiles, and the fraction of PAHs emitting the 
template is chosen to fit Spitzer IRS observations of local galaxies 
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II 



a 
Pi, 

3 



2000- 



1500 



1000- 



500 



q E os=0-75 
q E os=0-25 
fiducial [DOG 10] 
Smooth ISM 
no AGN 




0.6 0.8 
Time (Gyr) 



Figure 1. Evolution of 24/R ratio for our fiducial model (DOG10) while 
varying the ISM specification and inclusion of AGN. See text in ij |2.2.2l for 
details. 



jDopita et al.ll2005l ; IGroves et al.ll2008l: Ijonsson et alj|20ld) . Their 
absorption is calculated utilising the lli & Draind ( 120011) cross sec- 
tions. Varying this fraction changes the PAH luminosity (which can 
affect the 24 |im flux density at 2=2). We have performed test simu- 
lations exploring the range of PAH fractions, from to 1. Generally, 
the trends presented in this paper are insensitive to the inclusion of 
PAHs, though the normalisation of some trends can be. That is to 
say, the magnitude of the 24/R ratio is dependent on the fraction 
of emitting PAHs, though the evolutionary trends in the 24/R ra- 
tio are not. Tests investigating the limits (e.g. either 0% or 100% 
of carbonaceous grains with sizes a < 100A emit as PAHs) show 
that the observed 24 \im flux density of our model galaxies varies 
at most by a factor of ~2 (when observed at 2=2.0). We note that 
PAHs are not explicitly destroyed by the AGN. 

The evolving dust mass is calculated assuming a constant dust 
to metals ratio of 0.4, consistent with models and observations of 
both local U n iverse galaxie s , as well as those at higher redshifts 
dDweklll998l : I Vladiloi fl998l ; ICalura et alj|2008l) . The distribution 
and mass of interstellar dust is determined by the galax y evolution 
simulations. We use the IWeingartner & Draind J2001 ) dust grain 
model with R = Av /Eb-v = 3.15, as updated bv lDraine & Lil 
d2007h . This model include silicates and graphites. We note that 
silicate absorption features are not put in by hand into this model 
- rather their occurrence is the result of modeled radiative transfer 
effects. 

All photometric quantities in this paper are convolved with 
the appropriate instrument filter sensitivity function. Throughout 
this paper, all flux densities and flux density ratios come from both 
galaxies in the simulation box. While the simulation box is 100 kpc 
on a side, the nuclei of the galaxies (where the bulk of the emission 
arises from) never stray more than ~50 kpc from each other. This 
is roughly comparable to the Spitzer 24 \im beam at 2=2. Hence, 
during all the relevant stages of the interaction, both galaxies would 
appear as a single (blended) source in Spitzer 24pm surveys. 



2.2.2 ISM Specification and Other Free Parameters 

The ISM in our simulations is unresolved below the hydrodynamic 
resolution limit of ~ 100 pc. In order to model obscuration below 
these scales, we employ subgrid methods. In the remainder of this 
subsection, we discuss the effects of the subgrid assumptions on 
our final results. For the reader unfamiliar with the main results 
of this paper, it may be worthwhile returning to this section after 
reading through the main results to better understand the implica- 
tions of our subgrid techniques. In this section, we will refer to 
Figure [T] frequently, which summarises our comparison of varying 
ISM treatments. In this figure, we vary assumptions regarding the 
ISM and AGN for our fiducial merger, DOG 10 (with 80% initial 
gas fraction), and plot the evolution of the 24/R ratio. The thick- 
ness of each lightcurve owes to sightline variations in the SED. 

The young stellar clusters (with ages < 10 Myr) are as- 
sumed to reside in their nascent birthclouds, and have their spec- 
trum reddened accordingly. The birthclouds are comprised of HII 
regions and photodissociation regions (PDRs), whose sizes and 
temperatures are ca lculated with the 1 d imensional photoionisation 
code MAPPINGSIII ^Groves et alj|2008l). The HII regio ns evolve as 
one-dimensional mass-loss bubbles Icastor et"al]|l975l) and absorb 
much of the ionising UV radiation from the stellar clusters. In fact, 
the HII regions dominate the 24 \im continuum (observed-frame at 
2=2) flux density during starbursts. 

The time-averaged covering fraction of PDRs is a free pa- 
rameter, and can be related to the PDR clearing time scale as 
/pdr = exp{— t/tciear)- Here, we take the PDR covering fraction 
to be unity (/pdr = 1) which essentially means that massive O and 
B stars are assumed to be surrounded by their nascent birthclouds 
for the entirety of their lives (as the PDR clearing tim escale is then 
longe r than the lives of these stars; see Figure 6 of iGroves et al.l 
2008). While /pdr is an unconstrained parameter, some evidence 
exists that in local mergers a molecular interstellar medium with a 
large volume filling fac tor may blanket nuclear O and B stars for the 
major ity of their lives toownes & Solomon! 1998; Sakamo to et al.l 
1999). This suggests that a large PDR covering fraction may be 
a reasonable assumption in massive starbursts. The assumption of 
/pdr = 1 is conservative: removing the cold birthclouds typically 
increases the 24 |j.m flux density (observed at 2=2.0) (as emission 
from the HII regions can escape more unhindered). 

One potential issue is that the MAPPINGSIII SEDs are param- 
eterised in terms of their compactness, which is dependent both on 
the stell ar cluster masses an d the ISM pressure (see e.g. equations 
9-13 of IGroves e t al. 2008). The pressures seen in the central re- 
gions of the most massive high-2 mergers exceed the largest ones 
available in the MAPPINGSIII lookup tables by a factor of ~ 10- 100. 
This means that in reality the MAPPINGSIII HII/PDR SED will be 
hot ter than those at th e saturated values used here (see, Figure 5 
of IGroves et aUl2008l) . Because the stellar/HII/PDR SEDs domi- 
nate primarily during the peak starburst phase of a merger, the net 
effect of this is that we are likely underestimating the duty cycle 
of the DOG phase during merger-driven starbursts in our simulated 
galaxy sample. Because our inclusion of the MAPPINGSIII SEDs 
errs on the side of more conservative models, we proceed with this 
caveat in mind. 

Our simulations have no information about the structure of 
the ISM below the resolution limit outside of a stellar cluster's 
birthcloud. Effectively, what this means is that the radiative trans- 
fer models here do not account for absorption by GMCs outside of 
a given stellar cluster's birthcloud ( Jonsson et al. 2010). To place 
an upper limit on the potential attenuation from the cold ISM, we 
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have run test cases where we assume that each grid cell is uniformly 
filled with dust and gas with the mass returned from the SPH sim- 
ulations. This effectively assumes that all dust that would be as- 
sociated with GMCs in the ISM has a volume filling fraction of 
unity. Because the mass in the birthclouds is 'borrowed' from this 
gas mass budget, in this special scenario we require that /pdr = 
in order to avoid potential double-counting issues. In this extreme 
case, the ISM is optically thick at mid-infrared (rest-frame) wave- 
lengths, and the SED has a larger proportion of its flux at longer 
wavelengths. The purple curve in Figure [T] shows the evolution of 
the relatively depressed 24/ R ratio for this "Smooth ISM" case. We 
note that while the 24/R ratio is in this Smooth ISM model is not 
as large as the other cases, the galaxy's SED is still quite red. In 
this case, a larger fraction of the power in the SED is shifted into 
the FIR; as such, this galaxy would likely be selectable by e.g. Her- 
schel surveys. 

Our choice of equation of state may have an effect on the star 
formation history and structure of the IS M in the model galax- 
ies as well. Nominally, we utilise the full ISpringel & Hernquistl 
d2003h multi-phase EOS which simulates the pressurisation of the 
ISM owing to supernovae heating. This is somewhat stiffer than 
an isother mal EOS (the stiffest E OS is shown by the solid curve in 
Figure 4 of ISpringel et ai]20 05b, the isothermal EOS is the dashed 
line). Referring to this figure, the effective EOS employed can be 
quantified by a 'softening parameter' such that gsos = 1 repre- 
sents the stiffest EOS which aims to account for SNe pressurisation 
of the ISM. Similarly, gEOS = represents an isothermal EOS. 

Because we use the stiff EOS (gsos = 1), the gas is relatively 
pressurised which retards fragmentation. This affects both the star 
formation history as well as the structure of the ISM (galaxies with 
a softer EOS will be more clumpy which may affect optical depths). 
To investigate the consequences of our choice of EOS, we have 
examined the evolution of a fiducial galaxy merger, DOG 10 with 
softer EOSs (gEOS = 0.25 and gEOS = 0.75, linear interpolations 
between isothermal and gEOS = 1). 

In Figure Q] we show the 24/R lightcurves for our fiducial 
model, with the aforementioned modifications to soften the EOS 
(black and blue curves). The general effect of the softening of the 
EOS is to cause the ISM to become more unstable to fragmenta- 
tion, and reduce optical depths toward the central regions. Because 
of this, the observed SED is hotter for the gEOS = 0.75 and 0.25 
cases than the fiducial DOG 10, and the peak 24 [im flux density a 
factor of a few higher. Softening the equation of state increases the 
likelihood that a given simulation will be selectable as a DOG. 

There is some evidence that discs at hi gh redshift may 
indeed be clumpy. Both observations (e.g. Elmegreen et al. 
2007l;lDaddi et alj2008l.lGenzel et alj2006l ; lElmegreen et alj2009t 
Forster Schreiber et al. 200 9]) and models (e.g. iBournaud et al.l 
20081: iDekel et alj 120091: lElmegreerj 120091 : ICeverino et ail l20ld : 
Elmegreen & BurkertH2010T) suggest that disc galaxies at high-z 
contain relatively large clumps of gas. However, despite this, we 
choose to utilise the simulations with the stiffest equation of state 
(<Zeos = 1)- The reasons follow. First, with respect to simulat- 
ing the formation of DOGs, this is the middle-of-the-road choice 
amongst our ISM assumptions. For example, examination of Fig- 
ure Q] shows that while a softer EOS more easily produces DOGs, 
a smooth ISM supresses them. The model with birthclouds (and 
no other cold GMCs) and a stiff EOS has a test lightcurve which 
falls in between the two extremes. Second, we can constrain the 
models somewhat by the gas fractions at the time of merger. If we 
assume that mergers represent a reasonable mo del for the formation 
of high-redshift Submillimetre Galaxies (e.g. iTacconi et alj|2006l 



l2008tlNaravanan et alj2009ll20ld though see Dave et al. 2009 for 
counter-arguments), then the gas fractions in the model galaxies at 
the time of merger/peak submm flux density sho uld be compara- 
ble to those inferred from observations. Indeed, [Narayanan et al.l 
(2010) showed that the gas fractions at peak coalesc ence were in- 
deed ~ 4 0%, comparab l e to t he measurements of iBouche et al.l 
(2007) and lTacconi et al. | j2008l) . However, because our simulations 
do not include the infall of primordial gas from the IGM, there is 
no way to replenish the gas once it is consumed by star forma- 
tion. Thus, by utilising a softer EOS and increasing the level of disc 
fragmentation/star formation, the models quickly run out of gas and 
render the galaxy relatively gas poor during the peak submm phase. 
Tests show that the simulation with gEOS = 0.75(0.25) have gas 
fractions of ~ 20(10)% during their final coalescence/peak-submm 
phases. In this sense, our usage of a stiff EOS is a requirement 
to match certain physical properties of the most extreme high-z 
ULIRGs, necessitated by limitations of the models. Moreover, in 
the interest of making a more direct comparison with our previous 
attempts at modeling high-z ULIRGs (to wit, the SMG population 
iNaravanan et alj2009ll2oTol) , we utilise the stiffest EOS and birth- 
cloud model to remain consistent with the initial conditions chosen 
for those models. 

To summarise our exploration of the effects of ISM substruc- 
ture on the radiative transfer, we expect dumpiness of the ISM to 
exist on three scales: on the scale of birthclouds around young stel- 
lar clusters; in other GMCs outside of the birthclo uds; and on ^kpc 
scales as se en in observed discs at high-z (e.g. lElmegreen et al.l 
12001 2009). Our inclusion of the MAPPINGSIII calculations ac- 
counts for some obscuration on the birthcloud scale, and our pa- 
rameter search with the EOS has tested the possible effects of larger 
clumps on the radiative transfer. With respect to obscuration by 
GMCs outside of the parent birthcloud, the current generation of 
simulations only allow us to test the limiting cases of a smooth ob- 
scuring ISM with a volume filling fraction of unity. The results of 
these tests on the models are summarised in FigureQ] 

Finally, we can explore the role AGN play in our model for 
DOG formation. AGN heat the surrounding dust and contribute to 
rest frame 8(J.m thermal emission. This both increases the observed 
24/R ratio, as well as the duty cycle over which the galaxy has 
extremely red colors. The effect of not including AGN can be seen 
in the red curve in Figure Q] 



3 FORMATION AND EVOLUTION OF DOGS 
3.1 Physical Evolution of a Merger 

The formation of ULIRGs through mergers ha s a long history (e.g., 
ISanders et all 19881: ISanders & Mirabel 19961) . We therefore begin 
our investigation of DOGs with the ansatz that DOGs originate in 
gas rich galaxy mergers at high redshift. We will then go on to 
explore under what circumstances DOGs may be represented by 
non-merging galaxies (§ |4), and show that as observations probe 
sufficiently lower 24 \im flux density limits, less extreme condi- 
tions are sufficient to form DOGs. 

We first briefly orient the reader as to the generic physical 
evolution of a high-redshift binary gas rich galaxy merger. Here, 
we focus on fiducial model DOG 10 (Table [Q. DOG 10 is a 1:1 
Mdm « 6 x 10 12 M merger. Figure [2] shows the evolution of 
the morphology of DOG10 as viewed from a single SUNRISE cam- 
era. The boxes are 100 kpc/h on a side, and the galaxy is observed 
at z=2 (with the observed 24 [im, R and _B-band emission mapped 



© 2009 RAS, MNRAS 000,[T|-?? 



The Nature of z ~ 2 Dust-Obscured Galaxies 7 




Figure 2. Morphology of fiducial merger, DOG 10, modeled at redshift z=2. Colors are (observed-frame) MIPS 24 |xm, i?-band, and B-band mapped onto red, 
green and blue. Each panel is 100 kpc/h on a side, and the time stamp, 24 \aa flux density, and 2AIR ratio is denoted in each panel. The color of the 24/i? label 
represents if the mid-IR SED exhibits a bump or power-law nature (to be discussed in § 13. 3h where blue labels show bump galaxies, and red labels PL galaxies. 
Mergers naturally form DOGs. As the merger progresses, the galaxy becomes redder owing to increased 24 \aa flux density powered by the dust-enshrouded 
starburst. Upon final coalescence, when the SFR peaks, the galaxy is red enough that it may be observed as a DOG (24/i?> 1000). Subsequent contribution 
from a growing AGN will allow the galaxy to transition from a bump DOG to a PL DOG. 
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onto red, green and blue). DOG 1 is a very similar model to as the 
fiducial SMG (SMG10) studied in lNaravanan et al.l d2009h . We will 
present a more detailed comparison between DOGs and SMGs in 

In Figure [3] we show the evolution of the SFR, black hole 
accretion rate, and bolometric luminosity of DOG 10. The first pas- 
sage of the galaxies induces a ~200 MQyr -1 starburst. This ele- 
vated star formation rate is largely sustained for ~5x 10 8 yr while 
the galaxies inspiral toward final coalescence. As such, by the time 
of coalescence, the ga laxy system builds up a ^lO 11 Mq bulge 
dNaravanan et al.|[2010h . 

When the galaxies approach for final coalescence, tidal 
torquing on the gas funnels large quantities i nto t h e nu- 
cleus of the merged system (Barnes & Hernquist I l99l[ 1 19961 ; 



cleus of the merged sy stem (Barnes & Hernquist 1991, 1996; 
iMihos & Hernquistlll994 1 19961 ; ISpringel et al.N2005bl) , inducing 



a massive ~1000 MQyr~ starburst. Concomitantly, inflows fuel 
central black hole accretion. In the mergers modeled here, the black 
hole accretion rate can approach ~l-2 MQyr -1 . The associated 
AGN feedback contributes (along with gas consumption by star 
formation) to terminating the starburst, and may render the black 
hole optically luminous for a short period of time as feedback and 
gas de pletion clear sightlines (Springel et aL 2005b; Hopki nTet al.l 
2005 a b), a ccounting for th e observed population of bright quasars 
(Hopkins et al. 2006, 2008). The contribution of newly formed stars 
and an embedded AGN drive the ~10 13 Lq bolometric luminos- 
ity during the merger. It is during the final merger, when the star- 
burst, black hole accretion rate and bolometric luminosity are all 
near their peak (T ~ 0.65-0.75 Gyr; i.e., the shaded region in Fig- 
ure [3} t hat we expect the infrared luminosity of t he model galaxy 
to pea k djonsson et alj 2006; Yo unger et alJl200^ ; lNaravanan et al.l 
[20l3). 



3.2 The Evolution of the 24/ R Ratio in Mergers 

In Figure [4] we plot the evolution of the 24/R ratio for fiducial 
merger DOG10 placed at 2=2. As a reminder, Figures [2] and [3] 
serve as a reference for the morphology, black hole accretion rate, 
SFR, and bolometric luminosity for this galaxy. The grey shaded 
region shows the range of 24/R ratios as seen from 8 isotropically 
placed cameras outside the model galaxy, and thus reflects the de- 
pendence of the 24/R ratio on viewing angle. The yellow shaded 
region shows when the galaxy could be detected above the nomi- 
nal i*24>300 /J,]y detection threshold at z=2. We plot the 850 \im 
lightcurve as well (with units on the right ordinate), though defer 
discussion of this aspect of the figure until § 15. It . We additionally 
note when the galaxy would be a bump galaxy, and when it would 
be a PL galaxy (as the red/blue colors in the 24/R lightcurve), but 
defer discussion of this aspect of the galaxy's evolution to § 13.31 

Massive mergers at high redshift naturally produce DOGs dur- 
ing their final coalescence. The galaxy appears the reddest during 
the starburst/ AGN activity associated with final coalescence (in our 
fiducial DOG, we note that it is formally above the 24/R> 960 
threshold for only a brief period of time. More massive models [e.g. 
DOG1] form DOGs for longer duty cycles. DOG10 is, to some 
degree, the rough minimum galaxy mass necessary in our simu- 
lations to form a DOG). The observed 24 urn thermal emission 
comes from two sources. During the inspiral phase, warm dust in 
the diffuse ISM dominates the 24 \im emission. At times of elevated 
SFR (e.g., during final coalescence; T ~ 0.65 Gyr), HII regions ad- 
ditionally contribute a substantial fraction of the observed 24 urn 
thermal emission allowing the galaxy to be selected in 24 |j.m sur- 
veys. As gas consumption and the wind from the AGN terminate 
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Figure 3. Evolution of star formation rate, black hole accretion rate, and 
bolometric luminosity for fiducial merger model DOG 10. The bolometric 
luminosity includes contributions from both stars and black holes. The yel- 
low shaded region shows when the galaxy would be selected in a F24 > 300 
^iJy flux limited survey (the gap in the yellow region denotes a short time 
this particular simulated merger falls slightly below the i*24> 300/iJy flux 
density limit). The SFR, black hole accretion and luminosity all peak near 
when the galaxy is undergoing final coalescence. It is during this time that 
the galaxy system is expected to peak in infrared luminosity. We note that 
the black hole accretion rate and SFR are averages over 10 Myr intervals. 



both the starburst dSpringeletalj r2005a) and black hole growth 
(Youn ger et alj|2008bh , the 24/R ratio fades. Soon thereafter, the 
24 prn flux density drops (with the bolometric luminosity of the 
galaxy) such that the galaxy would no longer be detectable with 
current instruments. 

The 24/R> 1000 DOG criterion selects star-forming galaxies 
which have large quantities of obscuring dust in the diffuse ISM. 
The copious rest-frame UV photons produced by the starburst are 
reprocessed into rest-frame 8 [im thermal radiation by dust both 
within the HII regions and in the diffuse ISM. Though this point 
may seem trivial, it is to be underscored that z « 2 star-forming 
galaxies without large columns of dust in the diffuse ISM may still 
be selected as 24 p.m-bright galaxies (owing to the reprocessing of 
some of their UV photons in HII regions), but will not be reddened 
enough in their observed 24/i? to be selected as DOGs. We demon- 
strate this explicitly in Figure [5] where we show the 24 p.m and 
i?-band lightcurves for our fiducial model (DOG 10) both including 
diffuse dust (the fiducial model) and neglecting it. For the model 
neglecting diffuse dust, we still include HII regions surrounding 
stellar clusters. It is clear that the HII regions are responsible for a 
substantial fraction of the observed 24 \xm flux density. However, 
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Figure 4. 24/R and 850 \im light curves for fiducial model DOG10, mod- 
eled at z=2. The blue/red line is the sightline-averaged 24/R ratio as a 
function of time, and the color in this line denotes when the galaxy may 
be visible as a bump galaxy (blue) or PL galaxy (red). (See § 13.31 for a 
discussion regarding the transition from bump to PL mid-IR SED.) The 
grey shaded region highlights the range of 24/R ratios seen over 8 isotropic 
viewing angles, and thus reflects the sightline dependence of the 24/R ratio. 
The yellow shaded region shows when the galaxy's 24 |j.m flux density is 
above a fiducial limit of 300 fi]y. The break in the yellow shaded region is a 
short phase when the galaxy drops just below 300 /jjy. The black filled star 
and circle near the top show when the star formation rate and black hole 
accretion rates peak, respectively. The 850 urn flux density is overlaid as 
the purple curve, and has its units on the right axis. As the galaxy passes 
through its starburst dominated phase, it may be visible as both a Submil- 
limetre Galaxy (with -Fgso > 5 mJy), as well as a bump galaxy. The galaxy 
eventually transitions to becoming a PL galaxy, and may be either starburst 
or AGN dominated during this time. 



the diffuse dust is important in absorbing the UV radiation which 
escapes the HII regions; without this attenuation, the excess radia- 
tion which escapes the galaxy causes the normalisation of the 24/R 
ratio to fall significantly below the DOG selection criteria. 



3.3 Evolution of Bump DOG to Power-law DOG 

As discussed in § [TJ (see also D08), the mid-IR spectra of DOGs 
come in two flavors: those characterised by a 'b ump' at ~5 yim , 
and those exhibiting a power-law rise in the SED dDev et al.l2008T) . 
Bump features in the mid-IR SED of z ~ 2 galaxies are thought 
to originate from a combination of thermal emission from lower 
mass stars at rest-frame 1.6 urn, and emission from massive stars 
which show a local minimum in their atmospheric opacity at this 
wave length dJohnlll988l : ISimpson & Eisenhardlll 19991 : iFarrah et al.1 
2008). Similarly, the PL shape in the mid-IR spectra of z ~ 2 galax- 
ies is canonically associated with the dominant presence of AGN 
continuum (though, as we will show, the PL SED can originate in 
starburst dominated galaxies). 

In a merger-driven model for DOG formation, there is a natu- 
ral evolution in the mid-IR SEDs of z ~ 2 DOGs from bump-like 
to PL-like as the merger progresses. For the purposes of our anal- 
ysis, we fix our model galaxies at z=2, and define a bump galaxy 



Figure 5. 24 urn and i?-band lightcurves for the fiducial model (DOG 10) 
both including diffuse dust (solid lines) and neglecting it (dotted lines). We 
remind the reader that all models in this paper include diffuse dust with 
the exception of the dotted lightcurves in this figure. The models without 
diffuse dust still include HII regions surrounding stellar clusters. It is clear 
that the HII regions are responsible for a large fraction of the 24 urn flux. 
The diffuse dust is important in absorbing the UV radiation which escapes 
the HII regions; without this attenuation, the excess radiation which escapes 
the galaxy causes it to appear optically bright, and thus not selectable as a 
DOG. 

utilising a simple peak-finding methotjf] such that S3 < S5 and 
Ss > Ss (where S3, S5, and Sg, correspond to the filter-convolved 
fluxes from IRAC channels 1,2 and 4, respectively). Similarly, we 
define a PL galaxy as one where S3 < S5 < Ss. Returning to 
Figure [4] we now see that when we utilise this criteria for classi- 
fying the bump and PL phase of the galaxy's evolution, there is a 
transition during the merger from a bump galaxy to a PL galaxy 
(where the bump phase is denoted by the blue segment of the line, 
and the PL phase by the red one). Comparing with the morphol- 
ogy (Figure [2j, it is evident that this transition occurs during final 
coalescence. 

We can understand the details of the origin of the bump to 
PL transition more explicitly by examining the evolution of the 
SEDs of DOGs during the final coalescence stage. In Figure|6] we 
present the sightline-averaged SED from fiducial galaxy DOG10 as 
it evolves through its final coalescence phase. The observed SED is 
shown by the black solid line in each panel. The blue dashed SED 
represents the input STARBURST99 seeQ and the red dotted curve 
the input AGN template SED. Both the stellar and AGN SEDs are 
unreddened by diffuse dust in the galaxy (though the stars are red- 
dened by HII regions and PDRs). We list the time stamp of each 
panel for comparison with the morphology (Figure |2), and global 
24/R light curve (Figure^, and quote the 24/R ratio in each panel. 
The 24/R ratio is color coded based on if the DOG is a bump galaxy 
(blue), or PL galaxy (red). 

3 We note that this defin ition of "bump" ver sus "powerlaw" mid-IR SEDs 
is similar to that used by |Farrahet"ai] 1200811 . and is simpler to implement 
than the power-law fitting method used by D08. 

4 More specifically, the blue dashed line is the STARBURST99 input SED 
after it has passed through the HII regions and PDRs surrounding stellar 
cluster utilising MAPPINGSIII. 
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Figure 6. Model SEDs for 9 snapshots of fiducial galaxy DOG10 as it evolves through its final coalescence DOG phase. The black line is the sightline- 
averaged observed SED (at z=2). The blue dashed line represents the stellar SED (after passing through HII regions and PDRs), and the red dotted line is the 
unattenuated AGN spectrum. The SEDs are all redshifted to 2=2, and the abscissa is the observed-frame wavelength. The 24/R ratio is listed along with the 
time for each snapshot, and the color of the 24/R ratio represents when the galaxy would be selected as a bump galaxy (blue) versus a PL galaxy (red). A 
green vertical line is placed at the rough location of the redshifted stellar bump in the snapshots which are bump galaxy. We note that we label the galaxy as a 
bump or PL galaxy regardless of whether it classifies as a DOG by the fiducial 24/R criteria. The galaxy's DOG phase begins as a bump DOG (e.g., T = 0.67 
Gyr). As the merger progresses, the galaxy will be selected as a PL DOG, either owing to dust attenuation of the stellar bump (even at times when the AGN 
contribution is not particularly strong; e.g., T = 0.70 Gyr), or to direct contribution from the AGN continuum (e.g., T = 0.74 Gyr). It is important to note 
that even some sources that may appear as bump galaxies by eye (e.g. T = 0.69 Gyr) may not be classified as such when convolving with IRAC filters. 
Finally, we note that while only a single snapshot in this Figure is formally a DOG, the duty cycle increases with increasing galaxy mass. The fiducial DOG 
represents a rough lower limit to the galaxy mass which forms DOGs. 



In the early stages of the DOG phase (near the beginning of 
final coalescence, e.g., T = 0.67 Gyr), the mid-IR SED shows the 
characteristic stellar bump associated with a star-formation domi- 
nated infrared luminosity. As the merger evolves, the mid-IR SED 
becomes best characterised by a power-law shape. PL galaxies can 
arise from either starburst-dominated or AGN-dominated DOGs. 
Starburst-dominated PL galaxies owe their origin to large columns 
of dust obscuration associated with the final merging. To see this, 
consider panel T = 0.70 Gyr in Figure [6] During the merger, 
the intrinsic stellar flux density at A < 5 |i.m can be extinguished 
by up to an order of magnitude. This can force the SED to look 
power-law-like, even while the AGN contribution at these wave- 



lengths is small. The galaxy then evolves into an AGN-dominated 
PL galaxy. Here, the AGN SED contributes substantially to the ob- 
served SED at A > 8 urn, again giving a power-law mid-IR shape 
(e.g., T = 0.74 Gyr). As is clear, then, PL mid-IR SEDs from 
z ~ 2 galaxies may represent either starburst dominated or AGN 
dominated galaxies. The relationship between the AGN contribu- 
tion to PL mid-IR SEDs in z ~ 2 24 |j.m sources will be explored in 
more detail in a forthcoming paper (D. Narayanan et al., in prep.). 

Finally, we note that in the post-merger stage, the galaxy 
exhibits bluer colors with a bump-like SED. These ULIRGs are 
reminiscent of those detected which are less extreme than DOGs 
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Figure 7. Evolution of sightline averaged 2AIR lightcurves for the most 
massive two of the model isolated galaxies in Table[T](iDOGl and iDOG2) 
are shown as blue lines. Both gas fractions are plotted for iDOGl and 
iDOG2. The range of colors owing to viewing angle of the reddest galaxy is 
shown by the grey-hatched region. While isolated galaxies typically are not 
red enough to satisfy the DOG criteria, the most massive, metal-rich ones 
certainly can when in edge-on viewing angles. 



in deep Spitzer survey s lYan et all d2005l . 120071) ; lLonsdale et al.l 
j2009h : lYanetltl]j2010h . 



4 THE PHYSICAL FORM OF DOGS : ARE MERGERS 
NECESSARY? 

Thus far, we have seen that gas rich galaxy mergers at high red- 
shift provide a plausible means for forming luminous DOGs. Fur- 
thermore, this mechanism for DOG formation suggests a picture 
for connecting bump DOGs and PL DOGs. This leads to a natural 
question: are mergers a necessary prerequisite for the formation of 
DOGs? 

We begin with Figure [7] where we show the evolution of the 
24/i? ratio for the most massive isolated galaxies in Table Q] We 
plot the lightcurves for both gas initial fractions modeled - 40% 
and 60%. The lightcurves are sightline averaged, and we show 
the sightline-averaged disperion for the reddest galaxy through the 
grey-shaded region. Generally, isolated galaxies are not quite red 
enough to satisfy the fiducial 24/R> 1000 criteria for DOG se- 
lection. This said, the most massive, dust-rich ones can represent 
DOGs. For example, the 40% gas fraction, Mbar ~ 5 x 10 11 
Mq model (which starts with nearly solar metallicity in the closed- 
box model) approaches DOG-like redness. That said, the most mas- 
sive, dust-rich isolated galaxies can be selected as DOGs when in 
an edge-on configuration. 

We can further investigate the relative role isolated discs and 
mergers play toward z ~ 2 ULIRG formation in Figure [8] where 
we plot the SFR versus the 24 pm flux density (and rest-frame 8 
pm luminosity) at z=2 for every snapshot in our simulation sam- 
ple that has a 24/ R ratio > 1000. We discriminate both between 
mergers and isolated galaxies, as well as whether the galaxy is a 
bump or PL DOG. The larger symbols are galaxies which would 
also be detectable as SMGs, though we defer discussion of this as- 
pect until § 15.11 Mergers populate the bulk of the luminous (i<24> 



300 /^Jy) DOG population in this diagram. These systems, which 
are extremely gas/dust rich, become red enough to be selected as 
DOGs during close passages and final coalescence (e.g., Figure[4]l. 
Mergers provide an efficient means for generating the starburst 
and/or AGN activity to heat the dust sufficiently that the galaxy 
is visible at 24 pm above 300 pJy. Nominally, we find a minimum 
SFR > 500 MQyr -1 necessary when AGN activity is negligible 
to power F%i> 300 /^Jy sources. Similarly, for DOGs powered 
by AGN, we find minimum AGN accretion rates of A/bh > 0.5 
MQyr -1 are required to power the brightest DOGs. 

At lower 24 pm flux density limits, however, mergers are not 
always required for the formation of DOGs. Massive discs can be 
red enough to satisfy the requisite infrared-optical color ratio for 
DOG selection, though are less bolome trie ally luminous. 

The most luminous DOGs in Figure [8] te nd to all be PL 
sources, consistent with observations ( Dey et al. 200^). The iso- 
lated galaxies exhibit lower 24 pm flux densities, and are always 
bump DOGs. At lower flux densities, the SFR is correlated with the 
observed 24 pm flux density (C. Hayward et al., in prep.). As the 
contribution of AGN becomes non negligible, though, this relation 
breaks down (e.g., for most luminous sources). While the trends in 
Figure[8]are robust, we caution against interpreting the absolute val- 
ues of the 24 pm flux density too literally. Rest frame mid-IR PAH 
features contribute to the observed F24 at z ~ 2. Modulating the 
assumed fraction of carbonaceous grains which emit the PAH tem- 
plate can shift the observed 24 pm flux density of a given galaxy. 
Changing the PAH fraction has the strongest relative effect on the 
lowest flux density galaxies, and at its limits (e.g., a PAH fraction 
of or 1), can shift the observed flux density by up to a factor of 
~50% (though the shift is less for more luminous sources). 

To summarise, our models suggest that the most luminous 
DOGs (i.e., those with i<24> 300/iJy at z ~ 2) likely result from 
mergers. In contrast, the less luminous DOGs can result from a 
wider range of (less extreme) conditions, perhaps even from secu- 
lar evolution. Within this framework, there are a few salient points 
regarding the physical nature of DOGs that appear to be general to 
all DOGs selected at all flux density limits in our models. 

First, in the absence of any merging activity, a DOG will never 
go through a PL phase. Some merging activity is necessary to pro- 
duce the dust obscuration of UV radiation and/or AGN activity 
which drive the PL shape to the observed mid-IR spectra. Merg- 
ers are required to funnel sufficient gas to the central regions to 
suppress (or overwhelm) the stellar bump (and/or fuel the AGN). 
As such, PL DOGs constitute a relatively narrow class of galax- 
ies, whereas bump DOGs are significantly more diverse as they are 
well represented by both mergers and non-mergers. This is seen 
explicitly in Figure [8] 

Second, while it is not shown, there is a rough minimum 
galaxy mass necessary for forming DOGs. Nominally, galaxies 
with halos A/bar << 5 x 10 11 Mq will not have the requisite gas 
supply to generate the starburst/AGN activity that drives the ob- 
served 24/R ratio. When mapping from baryonic masses to halo 
masses (e.g. Table QJ, the halo masses ar e also consistent wit h 
the clustering measurements published by Bro dwin et al. (2008). 
Moreover, the halo masses are consistent with the inferred halo 
masses of high- redshift galaxie s of a comp arable luminosity - 
that is, SMGs dBlain et alj|2004l). and quasars JCroom et al]|2005i 



iLidz et al.ll2006l : lHopkins et alj2 007 ; Shen et al .1120071) 
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Figure 8. SFR versus 24 |xm flux density for all snapshots of all model 
galaxies which qualify as DOGs (24/ R> 1000, but with no 24 p flux 
density limit). On the top axis, we show the rest-frame 8 )^m luminosity. 
The different symbols denote whether the galaxies are mergers or isolated 
galaxies. The symbols which are double sized qualify both as DOGs, as well 
as Fgso > 5mJy, though we defer discussion of the SMG-DOG overlap un- 
til S I5.1I There is a clear relation between the 24 \im flux density and SFR 
for the bump DOGs, though the contribution of AGN can muddy this for PL 
galaxies. Generally, the most luminous 24 \im sources require merging ac- 
tivity, and have PL mid-IR SEDs. As galaxies become less bright, they tend 
toward bump (star formation-dominated) galaxies, and are well described 
by non-mergers as well. The reader should note that the gap in between the 
mergers and isolated galaxies is artificial, and betrays the limited parameter 
space investigated in these models. This gap represents the transition from 
isolated galaxies dominating the bulk of the DOG population to mergers. 
A broader range of models ran at more varied masses, merger orbits and 
gas fractions would fill this space in. It is moreover important to note that 
the flux densities on the abscissa are subject not only to the uncertainties 
present in this model (§|5J, but also the assumed redshift of 2=2.0. As such, 
the normalisation of this trend should be taken with the appropriate caution. 



5 RELATIONSHIP OF DOGS TO OTHER 
HIGH-REDSHIFT GALAXIES 

5.1 Submillimetre Galaxies 

A number of clues suggest an intimate link between hi gh-redshift 
850 |i.m- selected SMGs, and the DOG population dDev et all 
2009). First, circumstantial evidence indicates physical similari- 
ties between SMGs and DOGs. SMGs are thought to form in 
massive starbursts which both observation al and theoretic al ev- 
idence imply may be merger-driven (e.g., IChapman et alj 



Tacconi et al. 2006, 2008; Baugh et al. 2005; Swinbank et al. 



2003 



2008 



Narayanan et al. 2009, 2010, though see Dave et al. (2009) for an 



alternative interpretation). Our model for DOG formation sug- 
gests that luminous DOGs, as well, are likely formed at the fi- 
nal stages of a merger (§ |4). Similarly, the inferred SFRs in 
SMGs and PL DOGs appear to both be extremely large (~ 
700 - 1500 M^vr-^ISwinbanketall^OO^dKovacs etalfeooel; 



Menendez-Delmestre et al. 2007; Valiante et al. 2007; Dey et al. 



20081 ; IBussmann et al.l 2009; Desai et al. 2009). DOGs and SMGs 



cluster in similar mass halos lBrodwin etani2008l) . and tenta- 
tive evidence suggests that may they contain somewhat simi- 
lar mass stellar bulges (e.g. [Swinbank e t"al] |2004| ; iBorvs etHll 



2009; Michaio wskTet alj 12009). Second, 



iNaravanan et alj d2010f) 's merger-driven model for the formation 
of SMGs utilised similar merger models employed in this study. 
This overlap in modes of galaxy formation (a merger-based mode 
for both SMG and DOG formation) suggests not only a physical 
connection between SMGs and DOGs, but possibly an evolution- 
ary one. Finally, an d perhaps most con vincingly, an analysis of the 
GOODS-N field bv lPopeetai] d2008l) has found a bona fide over- 
lap between the SMG and DOG populations: Namely, some ~30% 
of SMGs have 24/R colors consistent with DOG selection. 

To assess the relationship between DOGs and SMGs, we now 
highlight the purple 850 |rm flux density lightcurve overplotted 
in Figure [4] As the galaxies merge (T ~0.65 Gyr), radiation 
from the ~ 1000-2000 MQyr" 1 starburst is intercepted by the cold 
ISM, and re-emitted by cold dust as submillimetre-wave emis sion 
dGroves eta l. 2008; Jonsson ~et al.l201(il ; lNara vanan etalj20ld) . At 
the same time, diffuse dust in the galaxy heavily reddens the opti- 
cal light, rendering the galaxy visible as a DOG. The model DOG10 
would be classified as a bump galaxy at the exact same time as it 
would be visible as an SMG (Figure |4). The formal DOG phase 
follows after the SMG phase in this model, though we note that 
numerous models in Table [T] show an overlap of the two phases. 
After the SMG phase, the fiducial galaxy transitions from a bump 
DOG to a PL DOG. Because the SMG phase and DOG phase of 
the galaxy are both triggered by the merger, the phenomena are 
intimately related. 

While Figure [4] is instructive as to a general sequence of 
events, it should not be construed to mean that all SMGs would 
be selected as DOGs, much less bump DOGs. Merger orbit, view- 
ing angle, and mass can shift the exact timing and magnitude of 
events to make the overlap less straightforward. To show this, in 
Figure [9] we plot the 850 |J.m flux density from every snapshot in 
every model in Table [T] making two selection cuts: the galaxy has 
to either be classified as an SMG (fgso >5 mJy), or a DOG (24AR> 
1000). We additionally note when the galaxies would be visible as 
bump sources, or PL DOGs. There is a significant overlap between 
the galaxy populations. A few features are apparent. 

First, while there is overlap, many SMGs will not be se- 
lected as DOGs, and vice versa. Second, SMGs may be ei- 
ther bump or PL sources. The fact that some SMGs have bump 
mid-IR SEDs makes sense as the 850 \im flux which SMGs 
are selected on is powered by ab sorption of stellar light by 
cold clouds ^Narayanan et alj bold) , and are generally thought 
to have their bolome t ric luminosity powered by star formation 
dChapman et al.ll2004l; lAlexander et atll2005l: Ivaliante et alj|2007l: 
Pope et al. 2008; Younger et al .1 12008 a; Menendez-Del mestre et alj 
2009). That SMGs can also be PL DOGs may seem counterintu- 



itive given that the dominant power source in SMGs is thought to 
be star formation. Recall, however, that some PL sources derive 
their mid-IR SED shapes from extremely heavy dust extinction at 5 
|j.m, even while the AGN continuum does not dominate the mid-IR 
flux density. While the PL DOG phase is roughly contemporane- 
ous with the rise in AGN continuum, they are not always exactly 
coincident in time. 

SMGs will typically overlap with the more (24 \im) luminous 
DOGs. This owes to the fact that in o ur model, SMGs derive from 
massive mergers (see Figure 2 from INaravanan et alJfeplOh . The 
overlap of SMGs with high 24 p.m flux density DOGs is apparent 
from Figure [8] where we highlight the symbols with double size 
which are galaxies which are selectable both as DOGs, as well as 
F 850 > 5mly DOGs. 

Finally, we caution that Figure|9]should not be interpreted as a 
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Figure 9. 850 \im flux density versus 24/R ratio for all 24 u.m-selected 
galaxies in our model sample. Blue Mangles show when the galaxies would 
be selected as bump DOGs, and red diamonds as PL DOGs. Only galaxies 
which are either SMGs or DOGs are plotted. The vertical and horizontal 
lines denote the DOG 24/R selection criteria and the SMG -F850 selec- 
tion criteria. While there is significant overlap between the SMG (i^gso > 5 
mJy) and DOG population, numerous DOGs are not detectable in the sub- 
millimetre, and vice-versa. It is important to recognise that the region of 
space where 24/R < 1000 and i*850< 5 mJy (shown by the grey hatched 
region) is not populated by construction, as in this space galaxies are nei- 
ther DOGs nor SMGs. In reality, a large number of galaxies reside in this 
portion of Fssq-24/R, space. 



measure of the relative fractions of SMGs which are DOGs or vice 
versa. To construct those statistics will require a convolution of the 
SMG and DOG duty cycles with cosmological galaxy merger rates, 
a task deferred to a future study. 



5.2 BzK Galaxies 

The BzK selection criteria have been successfully used to find both 
star-formin g and passively e volving galaxies atredshifts 1.4 < z < 
2.5 (e.g., iDaddi et alj|2004r) . These samples overlap the DOG and 
SMG redshift ranges, and therefore provide an interesting compar- 
ison. 

In FigureQO] we plot all DOGs in our simula tion sample on t he 
BzK color-color plot. In order to compare with lPope et al.l |2008)'s 
analysis of GOODS-N DOGs (selected with F 2i > 100 /LtJy) in BzK 
space, we include all DOGs with i*24> 100 fiJy at z= 2. We sep- 
arate bump and PL DOGs, and overplot the rough criteria for star- 
fo rming galaxies ((z — K)ab — (B — z)ab > —0.2) devised 
by Daddi etHEl J2004h (such that galaxies which lie above the line 
are star-forming BzK galaxies). Finally, we note the location of all 
galaxies in our simulation sample (whether or not they are DOGs) 
in (z — K) — (B — z) space by dots. 

Most DOGs in ou r sample are select ed as BzK galaxies. This 
result is not surprising. IDaddi et al.l i200"3) find that the majority of 
BzK galaxies in the GOODS-N field are 2 4 ^m detected. In more 
direct comparison to these model results, iPope et al.l J2008h find 
that nearly all DOGs in the GOODS-N field are selected as BzK 
galaxies. The resu lts here are in good agreement with the dataset of 
Popeet al. (2008), spanning a similar range in (B — z) and [z — K) 
colors. 



Figure 10. DOGs on BzK color-color plot. Bump and PL DOGs are de- 
noted, and the selection criteria for star-forming galaxies ((z — K)as — 
(B — z)a.b > —0.2) is shown as a solid line (star-forming BzK galaxies lie 
above the line). We note the location of all galaxies in our simulation sam- 
ple (whether or not they are DOGs) in (z — K) — (B — z) space by dots. 
For the purposes of comparisons with observational data sets, DOGs in this 
plot are selected with i<24> 100 /iJy. Most model DOGs are selectable as 
BzK star-forming galaxies, consisten t with studies of the overlap between 
GOODS-N BzK galaxies and DOGs JPope et alj2003) . 



There are two important caveats. First, while most DOGs can 
be selected as BzK galaxies, this does not mean that all BzK galax- 
ies will be selected as DOGs. Observationally, IPope et al.[ ( 2008) 
find that only 12% of BzK galaxies are selectable as DOGs. Sec- 
ond, even though most DOGs fall into the star-forming region of 
the BzK color-color plot, they may not be detectable at optical 
wavelengths. Because the BzK criteria covers a broad range of star- 
forming galaxies, the DOG selection criteria selects a dusty subset 
of the BzK's. 



6 OBSERVATIONAL CONSEQUENCES OF THE MODEL 

In order to assess the validity of our models, we compare the sim- 
ulated and observed SEDs. We then make testable predictions of 
these models, and provide analysis tools via public SED templates. 
When comparing to observations, we focus primarily on the dy- 
namic range of fluxes and color ratios. Matching the observed dis- 
tributions (e.g., luminosity functions, color distributions) requires 
cosmological galaxy formation simulations, and will be presented 
in forthcoming work. 



6.1 Comparisons to Observations: Optical-mm Wave SED 

Observed SEDs of DOGs with spectroscopic redshifts have been 
compiled by Bussmann et al. (2009, in prep). In Figure QT| we 
present the model SEDs for four merger-driven DOGs which span 
the (mass) range of galaxies which form F 2 4> 300/xIy DOGs 
(DOG1, DOG4, DOG7 and DOG10) and sample the full mass ra- 
tios of simulated mergers. We note that this is a simple mass (and 
mass ratio) scaling in simulated mergers, but aside from this, the 
initialisation of each merger (including galaxy orbit) remains iden- 
tical. The SED is modeled for every snapshot while the galaxy is a 
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DOG (above 7*24 > 300 ^iJy), and is averaged over all camera an- 
gles. The black line represents the time and camera-averaged SED, 
and the blue shaded region denotes the la temporal dispersion for 
each model. The SED is modeled to be at 2=2.0. The red squares 
show the mean of all data in the Bussmann et al. sample for DOGs 
above F 24 > 300/xJy within the redshift range 1.75 < z < 2.25, 
and the error bars denote the propagated errors. In order to more ac- 
curately compare with the Bussmann et al. sample, we performed 
a Monte Carlo drawing of the galaxies in our simulation sample 
that forced our final sample to be comprised of 50% bump DOGs 
and 50% PL DOGs in the composite SED. This is comparable to 
the breakdown of mid-IR SED types in the observational dataset of 
Bussmann et al. We note, however, that simply including all DOG 
snapshots in our simulation sample makes little difference on the 
final composite SED. W e additionally include template SEDs for 
local galaxies Mrk 231 dArmus et alll2007l : ICharv & Elbazll200ll. 
R. Chary, private communication), and Arp 220 dRieke et al. 2009), 
each scaled to match the Bussmann et al. MIPS 24 [im point and 
redshifted to 2=2.0. 

Generally, there is a reasonable match between the observed 
and modeled SEDs. The correspondence between the simulated 
optical, NIR and FIR SEDs, and those observed in DOGs im- 
plies that the models presented here have a similar UV flux from 
newly formed stars, stellar mass, and dust obscuration column as 
DOGs in nature. There are relatively minor discrepancies in the 
observed-frame NIR and mid-IR. The Mrk 23 1 and Arp 220 tem- 
plates both have weaknesses in matching the observed data. The 
Arp 220 template overpredicts both the rest-frame optical-NIR and 
FIR flux densities. Similarly, the Mrk 231 SED underpredicts the 
rest-frame UV substantially. Given the relatively improved match 
between these models and the observed data points (as opposed to 
the local templates), we provide the model SEDs as templates in 

Figure[T2]shows where the simulated galaxies fall on an IRAC 
color-color diagram. Here, we have taken our simulated galaxies, 
and redshifted them between 2:=1.5 to 2.4 incrementally (in steps 
of Az=0.1), and plotted all galaxies with a 24/7? ratio > 1000, 
and F24 > 300/iJy. The dots show the galaxies which have F24 
> 300/xJy, but do not have the requisite 24/7i ratio to be se- 
lected as a DOG. The PL DOGs show redder colors than the bump 
DOGs, which in turn are redder (in 5.8/3.6) than the non-DOG 
24jim sources. The simulated colors of all models lie roughly in 
the sa me region as the observations jDesai etld1l2009l ; lLacv et all 
l2004h . but span a s maller dynamic ran ge. The solid lines enclose 
the region noted by lLacv et alj (1 2004) to typically contain AGN. 
Both observed DOGs (Des ai et alj|2009l) and these models tend to 
fall primarily in the lLacv et al.1 {2004) AGN-wedge. While indeed 
DOGs powered by AGN fall into the AGN-wedge, many starburst- 
dominated galaxies appear to as well. 



6.2 Testable Predictions 

Our modeling suggests that the bright DOGs (7*24 > 300/iJy at 
2=2) are the result of mergers and represent a phase of massive 
galaxy evolution following a sub-mm bright phase. If true, the mod- 
els can be used to predict observational tests of this scenario. We 
note that these tests do not apply to the fainter DOGs, which are a 
more diverse population. 



6.2. 1 The Physical Masses of DOGs and Their Relation to SMGs 

In Figure [T3] we plot the range of black hole, stellar, dark matter, 
and H2 masses for all galaxies in our simulated sample at times they 
may be a DOG, discriminating both bump and PL DOGs (7*24 > 
300 /iJy). This prediction is made for galaxies that would be se- 
lected as DOGs between a redshift range of 2=1.6-2.4. We addi- 
tionally show the same prediction for all galaxies in our simula- 
tions which would be selected as SMGs (with Fsso > 5 mjy). The 
black hole, stellar, and halo masses are extracted directly from the 
hydrodynamic simulations. We note that the halo masses are pre- 
dicted by construction as the mass of the simulated galaxies is set 
by the initial halo mass. While the SPH simulations do not explic- 
itly follow the H2 molecular gas, we infer its mass by considering 
all of the dense, star-formi ng gas to be neutral, and ex tracting the 
H2/HI fraction utilising the Blitz & Rosolowsky j2006l) interstellar 
pressure-based prescription 1 Narayanan et alj|2009l) . 

Because we are considering a 7 7 24> 300 fily flux density limit, 
all of the DOGs in this plot formed from mergers. Mergers produce 
the bulk (~ 90%) of their s tellar mass prior to the final merging 
dNaravanan et al.l2008l.l200l . T.J. Cox et al. 2009, submitted). Con- 
sequently, DOGs which form at the final coalescence of a merger 
will typically have already formed the bulk of their stellar mass. 
As such, the bump DOGs and PL DOGs in our simulations tend to 
have relativ ely large stellar masses, M* > 10 11 , compa rable to ob- 
servations dLonsdale et al J 12009: Buss mann et "all l2009h . Because 
similar types of galaxies form SMGs and bright (7*24 > 300/iJy) 
DOGs (e.g. massive mergers), there is substantial overlap between 
the physical properties of bright DOGs and SMGs. Generally, only 
the halo masses are well-constrained for DOGs. As such, the re- 
maining physical properties (stellar, BH, gas masses) shown in Fig- 
ure [13] serve as a prediction. 

The growth story for black holes is markedly different from 
that of the stars. Black holes in mergers undergo the bulk of 
their growth during th e final merging (see, e.g., Figure 5 from 
iNaravanan et al J [20081) . As such, because the PL phase tends to 
follow the bump phase in these high-redshift merger models, on 
average, the black holes in PL DOGs will be larger than those in 
bump DOGs. 

We can utilise this growth history of stellar and black hole 
masses in mergers to provide a test for our evolutionary model 
for mergers to transition from an SMG phase to a PL DOG phase 
(e.g., Figure |4j- Merger remnants, in our simulations, lie on the 
MeH-ATbuigc relation owing to the termination of the starburst 
event by AGN feedback, and self-regulation of black hole growth 
dDi Matteo et alJl2005l:lRobertson et alfeOOdlHopkins et alj|2007l: 
lYounger et al. 2008b; Li et al. 2007). Because the black holes grow 
substantially during the final coalescence of a merger, one might 
expect that the black holes in SMGs will be smaller, on average than 
those in PL DOGs if SMGs are to evolve into PL DOGs as these 
models suggest (Figure |4j- In contrast, their stellar masses will be 
roughly comparable as the stellar masses are already in place at 
final merging. 

To test this, we compare the location of all bright SMGs 
and bright (F24 > 1 mjy) DOGs in our simulation sample on the 
MeH-ATbuige relation (Figure[T4]l. We bin the SMGs into average 
(7^850 > 5 mjy) and bright (7<850> 10 mjy) bins. The extremely 
bright (7*850 > 10 mjy) SMGs constitute a cleaner test as they are 
always galaxies at the peak of their starburst, whereas lower flux 
density sources (e.g., Fg$o~ 5 mjy) SMGs can be either lower mass 
galaxies at the peak of thei r starburst, or more luminous galaxies 
both pre and post-starburst dNaravanan e"tai]|2009l l201oh . Never- 
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Figure 11. Simulated SEDs redshifted to z=2 for model galaxies. The SED was calculated for all snapshots that satisfied the DOG selection criteria (24/i?> 
1000, F24> 300 /ijy) for four model galaxies which sample the full mass and merger mass ratio range of our simulation sample (DOGs 1, 4, 7 and 10): The 
lcr dispersion amongst these snapshots is shown by the blue shaded region, and the mean by the black solid lines. The model SED is compared to the compiled 
observational data points (red) by Bussmann et al. (in prep) for all F<2A> 300/^Jy DOGs with spectroscopic redshifts 1.75 < z < 2.25. If error bars are not 
shown, it is because they are smaller than the symbol representing the data point. In order to best compare with the Bussmann et al. sample, we did a Monte 
Carlo sampling of our model SEDs that enforced that we had 50% bump DOGs and 50% PL DOGs in the sample, comparable to the relative breakdown of 
mid-IR SED types in the Bussmann et al. sample. The magenta dotted lines and green dashed lines are template SEDs for local galaxies Mrk 231 and Arp 220, 
scaled in flux density to match the observed Bussmann et al. 24 |jm data point, and scaled to z=2.0. The modeled SED and observed data generally correspond 
well, suggesting a reasonable match between the modeled stellar mass and flux and obscuring dust columns and temperature in the modeled galaxies and 
DOGs in nature. That said, the modeled mid-IR SED shows some disagreement with the observed IRAC points (especially at 8 |rm). Neither the Mrk 231 and 
Arp 220 templates provide good matches to the observed data. 



theless, we include lower flux density (Fg5o> 5 mJy) SMGs for 
the purposes of comparison with existing data sets. 

Figure[l4]shows a direct prediction for our model. Most SMGs 
will lie below the present-day Mbk-M^\ ec relation as their black 
holes have significant growth t hey will undergo throu ghout the re- 
mainder of the merger (e.g., lAlexander et alj 12008). They typi- 
cally have stellar masses ranging from 10 11 — 10 12 Mq, with the 
more common SMGs tending toward the lower end of this mass 
rang^fl Bright PL DOGs, which occur contemporaneously with or 
later than SMGs in a merger's evolution (e.g., Figure[4} will range 



5 Inde ed comparison to the compilation of SMGs by Alexan der et alj 
1 2008) (the blue cross) shows that the average observed SMG is more sim- 
ilar to our SMG models which form lower mass galaxies. 



from falling below the present-day MBH-Afb u i gc relation, to di- 
rectly on/above it. 



6.2.2 T 0ust ofDOGs 

Forthcoming facilities which probe both sides of the peak of the 
SED of DOGs (e.g., Herschel, ALMA) will place strong constraints 
on the effective dust temperature. These dust temperatures will 
serve as a test for the predicted SEDs of DOGs from our models. 

In Figure [15] we plot the predicted dust temperature of our 
model DOGs as a function of their FIR luminosity (40-1000 \im). 
The dust temperature is derived by converting the rest-frame wave- 
length of the SED peak to a temperature (using Wien's displace- 
ment law), and is calculated for every snapshot where 24/R> 1000 
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Figure 12. Mid-IR color-color plot for all model galaxies which would be 
selected above a 24 \im> 300 /xJy flux density cut. DOGs are shown as 
colored points, with bump DOGs as blue symbols and PL DOGs as red 
symbols. Black dots represent 24 |xm selected galaxies which do not qualify 
as DOGs. The model galaxies are incrementally redshifted from £=1.5-2.4, 
and plotted when they are selectable as DOGs. The solid lines enclose the 
typical regi on of color-color s pace occupied by AGN-dominated galaxies 
as found bv lLacv et alj i2004|). The locus of PL DOGs lie in the top-right 
region of these plots (e.g., lLacvetal.l l 2004) owing to their rising mid-IR 
SEDs. The mod el DOGs lie ro ughly in the same region of the mid-IR as 
observed DOGs iDesai et all2009l) . though span a smaller dynamic range. 
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Figure 14. Bright (F 2 a> 1 mJy) PL DOGs and SMGs (in two flux den- 
sity bins) on the present day MeH-^bulgc relation. The purple circles 
denote bright SMGs, the green diamonds the avera ge SMGs, the blue 
cross the results from observed SMGs by Alexa nder et alj j2008), the red 
asterixes bright DOGs, and the grey shaded region the locally observed 
A^BH-A^bulgc relation (including dispersion). In our model, SMGs typ- 
ically precede the PL DOG phase in mergers. During this time, while the 
stellar mass is roughly in place, the black holes undergo significant growth, 
placing SMGs and PL DOGs at different locations on the MeH-A^bulgc 
relation. Observations pinpointing the location of SMGs and PL DOGs on 
the MeH-A^bulgc relation will serve as a test for our model connecting 
SMGs and DOGs. 




Predicted Masses (M®) 

Figure 13. Predicted range of physical properties of DOGs and SMGs 
(-F24 > 5 mJy). Shown are black hole masses, H2 masses, stellar masses 
and halo masses for galaxies which would be selected as DOGs between 
redshifts 2=1.6 and 2.4. Blue arrows denote bump DOGs and red arrows, 
PL DOGs. The vertical lines are simply to guide the eye to the axes. DOGs 
are relatively massive galaxies, with PL DOGs tending toward the higher 
mass end of the mass distribution. It is important to note that owing to rel- 
ative abundances, most observed DOGs are likely to occupy the lower end 
of the predicted ranges. Because similar types of galaxy (mergers) repre- 
sent SMGs and bright DOGs, they share an overlap in physical properties. 
Because only halo masses are known for DOGs, the remaining physical 
properties serve as a strong testable prediction of these models. 



and i<24> 100 /iJy (we utilise the lower flux density limit for com- 
pleteness). For reference, we show all times when the model galax- 
ies would be visible as SMGs (Fsso> 5mJy as well). While de- 
riving Tdu S t from the full SED is a different approach than that 
typically used by observers (who fit the observed data with SEDs 
varying Td us t and the dust spectral index, f3, which are degener- 
ate with one another), we prefer this method: Utilising the entire 
model SED avoids any assumptions about /3, which is observed to 
fall within a relatively large range of j3 = 1 — 2 for z ~ 2 SMGs 
and DOGs dKovacs et al.ll2006l : lBussmann et alj|2009h . 

For star-forming galaxies (including bump galaxies and 
SMGs), there is a general trend for brighter sources to have 
hotter dust temperatures (which is largely a selection effect for 
galaxies selected for cold dust emission). PL DOGs typically 
have hotter dust temperatures than bump DOGs and SMGs. 
These features are co nsistent with observat i ons of both high 
redshift galaxies (e . g.. iBussmann et alj |2009l : ICasev et"aT] 1 20091 : 
lYounger et aHl2009l;lKovacs et alJbOOaboirj) and local ULIRGs 



(e.g.. lSanders&Mirabelll 996). For example, Fig ure |15 | also shows 
z ~ 2 bump ULIRGs from Youn ger" et alj ( [20090 and lKovacs et alj 



(2010), SMGs from Kovacs et al. ( 2006), and two lower limits from 
PL DOGs from lBussmann et alj ((2009) as black circles, purple tri- 
angles, and cyan dia monds, respectively. The ULIRGs from the 
lYounger et al. (2009) sample were selected based on their stellar 
bump llluang et al. 2009), and are th us bump galaxies. For the 
bump ULIRGs of iKovacs et alj (2010), we divided the luminos- 
ity by a factor of 2 to convert from Lir to L40-1000, based on the 
approximate conversion for Mrk 231. Both the bump ULIRGs and 
SMGs span a similar range in Td us t-LiR space as these simulations. 
Dust temperature determinations for bright PL DOGs will provide 
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a more stringent test of these models as they are expected to show a 
larger range of dust temperatures than either SMGs or bump galax- 
ies. 



6.2.3 The CO Line Widths from DOGs 

iNaravanan et alJd2008ll2009l) found that CO traces the global kine- 
matics of a galaxy. The CO line widths can constrain the masses 
and dynamical states of galaxies which form DOGs, and thus con- 
stitute a testable prediction of these models. For example, while in 
an isolated disc galaxy, the line width is reflective of the virial ve- 
locity of the galaxy. In a merger, the CO linewidths may increase 
by a factor y/2 — 2, owing to multiple galaxies being in the obser- 
vational beam and the temporary disruption of molecular discs. In 
this picture, galaxies which owe their origin to mergers will exhibit 
larger CO line widths than secularly evolving galaxies. 

In order to examine the predicted CO line widths from DOGs, 
we employ TURTLEBEACH, a 3D non-local ther modynamic equi- 
librium molecula r line radiative transfer code dNaravanan et al.l 
12006. 12008, 2009). TURTLEBEACH is an iterative Monte Carlo code 
which considers both radiative and collisional processes in deter- 
mining the molecular excitation (and, hence , source functions), and 
assum es statistical equilibrium in the levels dNaravanan et al ■I2008L 
2009). For galaxy-scale simulations, we implement a subgrid pro- 
cedure for including a mass spectrum of giant molecular clouds 
as singula r isothermal spher es, which follow the Galactic mass 
spectrum (Rosolowskv 2005) and mass-radius relationship (e.g., 
ISolomon et alJI 1987b . Once the level populations have converged 
from the Monte Carlo code, the equation of radiative transfer is 
integrated through the grid to build the emergent spectrum. 

In Figure [H we plot the range of CO (J=3-2) FWHM for 
all model galaxies in Table Q] when they classify as a luminous 
(F24> 300 fi]y) DOG at z=2. As a comparison, we plot the CO 
FWHM from the same model galaxies when they would classify 
as an SMG (fs50> 5 mjy). The line widths are viewing angle de- 
pendent, and, as such, we include the line widths from 100 ran- 
dom sightlines for each model galaxy. The ranges are plotted in 
histogram form, with each random sightline of each snapshot qual- 
ifying as an individual galaxy for the histogram. It is crucial to 
note that the histogram does not represent a probability distribution 
function. Rather, the ordinate values are meaningless as the galax- 
ies are not drawn from a cosmological simulation representing true 
mass functions and merger rates. The robust quantity to take from 
Figure[T6]is the range of predicted CO line widths for DOGs and 
SMGs. 

The line widths from DOGs show a large range, from -~100 
km s _1 to > 1500 km s - . This owes to the fact that both rela- 
tively minor mergers (e.g., DOG9) as well as massive major merg- 
ers (e.g., DOG1) form DOGs. The former are in a relatively dy- 
namically relaxed state, and can be viewed along some viewing 
angles in a face-on configuration, causing small line widths. The 
latter almost always have relatively large line widths. The CO line 
widths from SMGs, in contrast, are typically quite large, ranging 
from ^500- 1500 km s" 1 . This owe s to their origin in massive ma- 
jor mergers (Naravanan et al. 2010). Hence, DOGs which may be 
in a disc-like configuration at times, are expected to espouse a much 
larger range of CO line widths than SMGs. Tentative o bservational 
evidence suggests that this may indeed be the case dGreve et"al] 
l2005l ; lTacconi et alj200dl2008l ; lYan et alfcoid) . 
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Figure 16. Predicted CO (J=3-2) line widths (FWHM) for model galaxies 
when they qualify as DOGs and/or SMGs. The latter galaxy population is 
selected via a submillimetre flux density cut, Fgao > 5 mjy. Line widths in 
galaxies are a function of their mass and dynamical state. Because DOGs 
can be either minor mergers (which may maintain a disc-like gas configu- 
ration), or major mergers (which are dynamically hot), the CO line widths 
span a large range (FWHM ~ 100—1500 km s — 1 ). In contrast, because our 
modeled SMGs all o we their origin to massiv e mergers, they typically have 
large CO line widths (Narayanan et al. 200<J). (See text for further details.) 
It is important to note that the histograms are not true distributions, but 
rather simply the range of expected CO line widths from these galaxy pop- 
ulations. Because line width roughly scales with mass, the average galaxy 
will have line widths toward the lower end of the predicted ranges. 



6.3 Analysis Tools: IR SED Templates 

We provide public SED templates of z ~ 2 DOGs both to serve as 
a testable prediction for these models, as well as an analysis tool 
for forthcoming JWST, Herschel, 1.1mm and ALMA observations. 
The templates are at model spectral resolution covering 968 
wavelengths between approximately 0.3 micron and 3 mm. The 
SEDs are provided at z=2.0, the mean redshift for DOGs, and are 
calculated utilising every model in the simulation sample. The 
templates are provided as a function of total infrared luminosity 
(1-1000 prn). The templates as well as an IDL script to read 
them in (and redshift them to arbitrary redshifts) may be found at 
http://www.cfa.harvard.edu/~dnarayan/DOG_SED_Templates.html 



7 DISCUSSION 

7.1 The Role of DOGs in Galaxy Evolution 

The models presented here suggest that the bright DOGs are the 
result of mergers, whereas the fainter DOGs can have diverse ori- 
gins. At the high flux density end (F 2 4> 300 ^Jy at z « 2), DOGs 
represent the coalescence phase of massive, gas rich galaxy merg- 
ers residing in large (A/dm ~ 10 13 Mq) halos. Bolometrically, 
these galaxies are amongst the most luminous at z ~ 2, and form 
stars at rates of SFR > 500 MQyr -1 . Forming bulges with masses 
a few x 10 11 Mq during the course of their mergers, these galax- 
ies likely are precursors to massive present-day spheroids. These 
merger-driven DOGs share an intimate link with SMGs, and repre- 
sent objects actively transitioning between having their IR luminos- 
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Figure 15. Predicted Td us t for DOGs above i<24> 100 fi]y and SMGs with Fg5o> 5 mjy as a function of Ltr(40-1000 |xm). The simulated Td us t is 
calculated via the peak of the sightline averaged SED. There is a general trend for blighter DOGs to exhibit hott er dust temperatures. PL DOGs hav e hotte r 
dust temperatures th an bump DOGs, as well as SMGs. For comparison, we show observational data points from Young eT et al.l fe009l) ; IKovacs et all fcOlOl) . 
iKovacs et al. 12006) and Bus smann et alj j2009) which represent 1 mm-selected z ~ 2 bump ULIRGs, z ~ 2 SMGs, and PL DOGs (lower limits) respectively. 
There is a reasonable agreement with the high-zobservations and these simulations. 



ity dominated by reprocessed stellar radiation to being dominated 
by a dust-enshrouded AGN. These massive mergers provide a nat- 
ural avenue for supermassive black hole growth, fueling the growth 
of ~ 10 9 M black holes. 

The bright DOGs (i<24> 300/iJy) typically are undergoing 
either a starburst (here, roughly >500 JViQyr -1 ) or contain an ac- 
tively accreting AGN (Mbh > 0.5 MQyr -1 ). These phenom- 
ena are most easily achieved via mergers. As was seen in Fig- 
ure [8] 1:1 mergers are not necessary, and more "minor" mergers 
(here, 1:3) can form bright DOGs. Given the increasing relative 
number densities of minor mergers (e.g., Fakhouri&Mij l2008l : 
Hop kins et ai]|2009h . it may be that these sorts of mergers domi- 
nate bright DOG formation. It is interesting to note that the higher 
mass-ratio mergers (e.g., 1:10) are not dissimilar qualitatively from 
galaxies which are accreting massive amounts of co ld gas from 
their parent halo at z ~ 2 (e. g ., "Stream Fed Galaxies" ; 
2005 ; iBirnboim et alj|2007l : iDekel et al ] |2008t 120091 : 



2009). 



Keres etal 



Keres et al 



At lower (24 ^m) flux densities (F 2 4<100 /uJy), DOGs can 
result from either galaxy mergers or galaxies undergoing secular 
evolution, not in any active stage of heightened star formation or 
AGN activity. These galaxies tend to have lower SFRs (-~50-100 



MQyr -1 ), though are of comparable mass as their merging coun- 
terparts. In order to be selected as DOGs, these galaxies must be 
viewed in an edge-on configuration which maximizes the dust red- 
denning. These galaxies may have overlap with the high redshift 
BzK population, and (in these simulations) always have their in- 
frared luminosity dominated by star formation. DOGs form a con- 
tinuum, then, of galaxies, parameterised by their 24 (xm luminosity. 
As such, there is, of course, overlap between merging and non- 
merging DOG populations in this sequence. 

Quantifying the exact number statistics of types of galaxies 
(e.g galaxy mass, merger mass ratio) in a given luminosity bin re- 
quires convolving light curves with cosmological galaxy merger 
rates, a task that is outside the scope of this paper. That said, these 
model results do point us toward some understanding of the phys- 
ical form of DOGs at various luminosities. Figure [8] demonstrates 
that mergers dominate the high flux density end, whereas isolated 
disc galaxies contribute to the fainter population. Genetically, we 
can say that the lower flux density DOGs, commonly identified in 
the deep, narrow field surveys, are likely to be a different breed 
than higher flux density DOGs which are typically uncovered by 
the wide-area Spitzer surveys. 

The idea that DOGs have luminosity-dependent physical 
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Figure 17. Fractional AGN contribution to bolometric luminosity for DOGs 
(above F24> 100 fi]y) versus 24 |j.m flux density. Generally, more lumi- 
nous DOGs tend to have larger AGN contributions to their total luminosity. 



properties has already been seen in observational data sets. For ex- 
ample, there is some evidence that more luminous DOGs are more 
likely to harbor AGN. D08 noted that at increasing 24 pm flux den- 
sities, a larger fraction of galaxies are X-ray detected and have PL 
mid-IR SEDs. Similarly, a number of studies have shown that lumi- 
nous DOGs ( J_4> 500/iJy) may be Compton thick and contain ob- 
scured AGN iFiore et alj2008ll2009l ; lGeorgantopoulos et al.120091 ; 
ISacchi et al.l2009T) . In our simulations, we see that brighter sources 
have increasing contribution to the bolometric luminosity by AGN 
(Figure [771). While assessing whether or not these AGN are Comp- 
ton thick is outside the scope of this work, we can make some crude 
estimates along these lines. If we assume that the molecular ISM in 
our simulations is uniformly distributed in th e nuclear regions (as 
is tho ught to be the case for local mergers; iDownes & Solomon 
1998), and not terribly clumpy, then we find some sightlines to- 
wards the nucleus may be marginally Compton thick (A^h ~ 10 24 
cm -2 ). However, a full investigation assessing the sightline depen- 
dence of column densities as well as the assumptions regarding 
the cold-phase clumping will be necessary to fully address whether 
these simulated DOGs are Compton thick. 

Further evidence for l uminosity-depend e nt phy sical proper- 
ties in DOGs comes from Melbour ne et al.l J2009h . who found 
an increasing concentration and decreasing half-li ght radius with 
increa sing 24 pm flux density in DOGs. Finally, iBrodwin et al.l 
found that DOGs exhibit luminosity-dependent clustering, 
such that the most luminous DOGs in their sample {F2i> 600/iJy) 
tend to reside in ~ 10 13 Mq halos, while less luminous (F24, > 300 
/ijy) DOGs are associated with lower mass A/dm ~ 2 x 10 12 
Mq halos. All of these are generally consistent with our picture for 
more luminous DOGs to be represented by high-redshift mergers, 
and less luminous DOGs displaying less star formation and AGN 
activity (Figure \V7\. 



8 SUMMARY AND CONCLUSIONS 

Utilising a combination of polychromatic radiative transfer calcu- 
lations and hydrodynamic simulations of high-redshift galaxy evo- 



lution, we have formulated a physical model for the formation and 
evolution of z ~ 2 Dust-Obscured Galaxies. We have utilised the 
very similar models as those employed in previous studies aimed 
at investigating the formation and evolution of a similar class of 
2 ~ 2 ULIRGS: the Submillimetre Galaxy population. 

While the uncertainty in the models is dominated by small- 
scale physics which is unresolved by the simulations, our method- 
ology allows us to make the following general conclusions re- 
garding the origin of DOGs, and their relationship to other high- 2 
ULIRGs: 

• DOGs are a diverse class of galaxies, ranging from secu- 
larly evolving star-forming disc galaxies (forming stars at ~50-100 
Mgyr" 1 ) to extreme gas-rich galaxy mergers forming stars at > 
1000 M yr -1 . 

• The most luminous DOGs (i<24> 300 /ijy) are well repre- 
sented by mergers of massive (Mom ~ 5 x 10 12 — 10 13 Mq) 
galaxies. These galaxies are actively transitioning from being star- 
burst dominated to being AGN dominated. At decreasing 24 pm 
flux densities (-F_4< 100 /^Jy), DOGs may either be galaxy merg- 
ers, or secularly evolving disc galaxies with more modest (~50-100 
MQyr" 1 ) SFRs. 

• Luminous, merger-driven DOGs naturally transition from 
having a bump-like mid-IR SED to a power-law (PL) shape, as the 
dominant power source transitions from star formation to the cen- 
tral AGN. That said, there is an overlap period where star-formation 
dominated sources can appear as PL galaxies. 

• The most luminous DOGs assemble both significant stellar 
masses (M* ~ 10 Mq), as well as contribute toward the growth 
of supermassive (Mbh ~ 10 9 Mq) black holes. 

• Merger-driven DOGs overlap with the 2 ~ 2 Submillimetre 
Galaxy population. SMGs generally represent the earlier, starburst- 
dominated phase of the DOGs evolution. This can be tested via the 
location of SMGs and DOGs on the MBH-A-buige relation. 

In advance of Herschel, JWST and ALMA, we provide the 
following testable predictions for our model of DOG formation: 

• We quantify the expected range of black hole, H2, stellar, and 
halo masses for luminous (i<24> 300 fily) DOGs as well as SMGs. 

• We provide dust temperatures for DOGs selected at various 24 
pm flux density limits. 

• We detail the location of SMGs and bright (F24 > 1 mjy) 
DOGs on the Merr-A/buigc relation as a test for the modeled evo- 
lutionary scenario that SMGs evolve into PL DOGs. 

• We provide a prediction for the CO line widths from DOGs, 
and suggest that they will be of order the largest observed line 
widths at high redshift, comparable to 2 ~ 2 SMGs. 

Finally, we provide model SED templates for 2=2 
DOGs as a function of total infrared luminosity (Lir= 
1-1000 pm). These are provided publicly, along with 
an DDL script for easily reading in the templates at 
http://www.cfa.harvard.edu/~dnarayan/DOG_SED_Templates.html 
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